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ABSTRACT 

We present high resolution observations of an extremely metal-poor damped Lycv 
system, at z a b s = 2.3400972 in the spectrum of the QSO J0035— 0918, exhibiting 
an abundance pattern consistent with model predictions for the supernova yields of 
Population III stars. Specifically, this DLA has [Fe/H] ~ —3, shows a clear 'odd-even' 
effect, and is C-rich with [C/Fe] = +1.53, a factor of ~ 20 greater than reported in any 
other damped Lya system. In analogy to the carbon-enhanced metal-poor stars in the 
Galactic halo (with [C/Fe] > +1.0), this is the first reported case of a carbon-enhanced 
damped Lya system. We determine an upper limit to the mass of 12 C, A/( 12 C)^ 
200 Mq, which depends on the unknown gas density n(H); if n(H) > 1 cm~ 3 (which 
is quite likely for this DLA given its low velocity dispersion), then M( 12 C) ^ 2 Mq, 
consistent with pollution by only a few prior supernovae. We speculate that DLAs 
such as the one reported here may represent the 'missing link' between the yields of 
Pop III stars and their later incorporation in the class of carbon-enhanced metal-poor 
stars which show no enhancement of neutron-capture elements (CEMP-no stars). 
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1 INTRODUCTION 

Damped Lya systems (DLAs) are the neutral gas reser- 
voirs at high redshift that have, by definition, neutral hy- 
drogen column de nsities in excess of 10 20 - 3 atoms cm (see 
' I Wolfe et af]|2005l for a review). At these high column den- 
sities the gas is self-shielded (e.g. IVladilo et al.|[200ll ). re- 
sulting in a simple ionization structure which facilitates 
the derivation of element abundances. Moreover, the abun- 
dances thus derived are independent of the geometrical con- 
figuration and thermodynamical state of the gas, and of 
most other fac tors which com plicate the analysis of stel- 
lar spectra (e.g. lAsplundl 2005) . The largest uncertainties in 
DLA abundance studies are due to the effects of line satu- 
ration and dust depletion (for some elements), although the 
latter of these concerns is found to be minimal when the 
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meta l licity of the DLA is below ~ 1/100 Zq ( Pettini et al.l 
1 19971 : iProchaska fc Wolfdl2002l lAkerman et ail 



20051 ). 



In recent years, these very metal-poor (VMP) DLAs, 
with metallicities [Fe/fij]]< —2, have attracted an increas- 
ing amount of interest because of their potential for probing 
gas which may still bear the chemical imprint of the first 
few generations of stars to have formed in the Universe (e.g. 
lErni et al.| [2006). They are thus an extremely valuable com- 
plement, at high redshifts, to local studies of the oldest and 
most metal-poor stars in the Galactic halo. 

The VMP DLA regime was largely unexplored until 
very recently, when it became possible to identify candi- 
date metal-poor DLAs in Sloan Digital Sky Survey (SDSS) 
quasars, and then measure their chemical composition with 
high resolution follow-up spectroscopy. The first high spec- 
tral resolution (R > 30000, FW HM < 10 km s" 1 ) s ample of 
VMP DLAs was compiled by iPettini et ail l|200Sl ), whose 
main goal was to compare the relative abundances of C, N, 



1 We adopt the standard notation: [A/B] = log(7VA/A r B) — 
log(A r A/A r B)0 where b refers to the number of atoms of ele- 
ment A,B. 
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and O with the trend s observed in VMP h alo stars in our 
Galaxy. More recently, iPenprase etafl (|2010h have presented 
medium resolution (FWHM ~ 60 km s _1 ) spectroscopy of 
a sample of 27 VMP DLAs to explore the general proper- 
ties of DLAs in the VMP regime. However, as Penprase et 
al. acknowledge themselves, there are difficulties with ac- 
curately measuring column densities from medium (as op- 
posed to high) resolution data, given that most DLAs with 
metallicities less than f/fOO solar exhibit very low veloc- 
ity dispersions , with metal line wid t hs les s than 10 km s~ 
jLedoux et all 120061 : iMurphv et "all 120071 ; iProchaska et all 
120081 ). Under these circumstances, line saturation can easily 
be overlooked. 

Perhaps the most startling result from abundance stud- 
ies of VMP DLAs are the near -solar values of [C/O] at low 
metallicity l|Pettini et al . 2008), in line with the peculiar up- 
turn in the [C/O] abundanc e below [O/H] — 1.0 in Galac- 
tic halo stars discovered by 



Akerman et all ||2004) and later 



confirmed bv lFabbian et al.l (1200981 ). I Akerman et al.1 (120041 ) 

attributed this behaviour to an i ncreased C yield from ear- 
lier generations of massive stars. IPenprase et al.l (|2010l ) ex- 
tended this work and reported a number of DLAs with su- 
persolar [C/O] suggesting that this ratio continues to in- 
crease with decreasing [O/H]. A more recent compilation of 
high spectral resolution observations of VMP DLAs, how- 
ever, suggests that t he ratio in fact plateaus at [C/O] ~ —0.2 
jCooke et al.ll2010bl ). 

In contrast to the relatively new interest in VMP DLAs, 
studies of metal-poor Galactic halo stars have received on- 
going attention for mor e than two decades, most recently 
from the dedicated HK ([Beers. Preston, fc Shectmanlll992i ) 
and HES (|Christlieb et al.ll200ll ) surveys. A relevant result 
emerging from this work is that nearly one-quarter of all 
metal-poor stars with [Fe/H] < —2.0 ex hibit a marked car- 
bon enhancement, with [C /Fe] > +1.0 (|Beers fc Christliebl 
120051 ; iLucatello et al. 2006). These are collectively known as 
carbon-enhanced metal-poor stars (CEMP stars), and have 
been subdivided into four classes based on the abundances 
of their neutron-capture elements: (i & ii) The CEMP-s and 
CEMP-r classes, with enhancements of elements produced 
predominantly by the s-process and r-process respectively; 
(iii) the CEMP-rs class, with enhancements in both the s- 
and r-process elements; and (iv) the CEMP-no class, which 
exhibits no such enhancements. 

For further details of these classes, and the likely ori- 
gi ns of their carbon enh ancements, we direct the reader 
to iMasseron et al.l (|2010h . In short, there is reasonable ev- 
idence to suggest that CEMP-s and CEMP-rs stars arc 
extrinsically polluted by a now extinct asymptotic giant 
branch (AGB) companion. The origin of the CEMP-no 
class, however, is not yet firmly established. Whilst sev- 
eral models invoking mass transfer from an AGB com- 
pani on could explain the lack of neutron-capture elements 
(e.g. Fuiimoto. Ikeda. fc I bcn 2000; Sies s. Gorielv. fc Langerl 



12004 ). radial- velocity measurements have not yet confirmed 
whether CEMP-no stars are host to binary companions. In- 
deed, the apparent difference in the metallicity distributions 
between CEMP-no stars and the other CEMP classes, in 
the sense that CEMP-no stars are more numerous at lower 
metallicity, might suggest that a mechanism other than the 
AGB binary transfe r scenario produc es the C enhancement 
in CEMP-no stars |Aoki et al.ll2007f ). In all likelihood, as 



pointed out bv lMasseron et all l|2010l ). there is a continuous 
link that connects some CEMP-no stars with CEMP-s stars, 
while the carbon enhancement of other CEMP-no stars may 
have a different origin. 

Models of core-collapse supernova yields from Popula- 
tion III (or near m etal-free) stars do enterta in high C yields 
relative to Fe (e.g. IWooslev fc Weaver|[T995l ). Population III 
enrichment is a particularly intriguing explanation for the 
origin of the carbon enhancement in some CEMP stars, since 
the three most Fe-poor halo stars all exhibit carbon enhance- 
ments relative to iron (|Christlieb et al l |2002| ; iFrebel et ajj 
120051 ; iNorris et al"1l2007l) . Moreover, the fraction of metal- 
poor stars that exhibit a c arbon-enhancement incr eases 
with decreasing metallicity (Be ers fc Christliebl 120051 ). In 
this Population III enrichment scenario, the carbon enhance- 
ment in the extremely metal-poor regime reflects the initial 
composition of the gas from which these stars formed, rather 
than resulting from mass transfer from an evolved compan- 
ion. 

Whilst the physics behind core-collapse supernovae 
is poorly constrained, several parameterised models have 
been developed to calculate the expected yields from zero- 
metallicity progenitors to compare with the observations of 
the most Fe-poor CEMP stars, as well as the more 'nor- 
mal' non-CEMP stars. The most important (and largely 
unknown) quantities are the degree of mixing during the 
supernova explosion and the a mount of fallback onto th e 
remnant black ho l e ther eafter l|Umeda fc Nomotol 120031 ). 
lUmeda fc Nomotol (|2003l ) parameterised both quantaties, 
then suitably selected the appropriate val ues that reproduce 
the ob served stellar abundance patterns. iHeger fc Wooslevl 
(2010), on the other hand, parameterise only the mixing pa- 
rameter by applying a running boxcar filter over the star 
following the explosion. This prescription successfully re- 
produces the chemical composition of the extr emely metal- 
poor ( non-CEMP) halo stars from the study bv lCavrel et al.l 
(|200J), as well as that of the most F e-poor CEMP stars 
HE0107-5240 (IChristlieb et al.1 12002|) and HE 1 327- 2326 
jFrebelet al.1 120051 ). IJoggerst, Wooslev. fc Hegerl (|2009l ) ex- 
tended this work by ma p ping the one-dimensional mod- 
els by IHeger fc Wooslevl l|2010l ) onto a two-dimensional 
grid to follow Rayleigh- Taylor induced mixing after explo- 
sive nu clear burning. Whilst the ljoggerst. Wooslev. fc Hegerl 
(2009) models are physically motivated, they are unable to 
reproduce the relatively high levels of nitrogen enrichment 
that are observed in the most Fe-poor stars, nor are they 
able to produce sufficient Fe-peak elements. 

These concerns a re alleviated with models that in- 
clude rotation (e.g. | Mevnet. Ekstrom. fc Maederl 120061 ; 
iHirschil [20071 ; iMevnet et al l |2010| ). since rotation induces 
additional mixing and mass loss in the pre-supernova 
phase of low- metallicity stars. The only simulation avail- 
able to date that investigates the effects of both ro- 
tational and Rayleigh- Taylor mixing, as well as incor- 
porating a realistic pr escription of fa llback, are those 
presented recently by Joggerst et al.l (|2010al ) (see also 
IJoggerst. Almgren. fc Wooslevll2010bl). By introducing rota - 
tion in their zero- metallicity models, [jj3gj2£rstTet al. (2010a) 
report an increased nitrogen yield, as well as an increased 
Fe-peak element yield, however, these models are not able to 
reproduce the high CNO to Fe ratios observed in the most 
metal-poor stars. In summary, it is still a matter of some 
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debate whether these extremely metal-poor early Popula- 
tion II stars: (i) were borne out of gas which had previously 
been carbon-enriched by the first stars; or (ii) received a 
CNO top-up from a companion star; and/or (iii) exhibit a 
degree of self-pollution from their own nucleosynthesis. 

The first possibility, that at least some CEMP stars 
were borne out of carbon-enhanced gas, is given additional 
support from the discovery reported here of an extremely 
metal-poor DLA at z a b s = 2.3400972 exhibiting a carbon en- 
hancement akin to that measured in Galactic CEMP stars. 
No other case was known until now, given (i) the rarity of 
DLAs with [Fe/H] < —2, whic h lie in the tail of the metallic- 
ity distribution of DLAs (e.g. iPontzen et ai1l2008h . and (ii) 
the difficulty in measuring [C/H] in DLAs, where the rele- 
vant absorption lines are often saturated even in the VMP 
regime. We speculate that this DLA may be the 'missing 
link' between the first few generations of stars and some of 
the CEMP stars in the Galactic halo. 

The paper is organized as follows. Section 2 summa- 
rizes the observations and data reduction. We analyse the 
absorption lines in the damped Lya system and deduce cor- 
responding element abundances in Section 3. In Section 4 
we scrutinize potential issues that could masquerade as a 
carbon enhancement. After ruling out these possibilities, we 
discuss in Section 5 possible origins of this enhancement, 
and consider the implications of our finding for models of 
CEMP stars. Finally, we summarize our results and draw 
our conclusions in Section 6. 



2 OBSERVATIONS AND DATA REDUCTION 

The m r = 18.89, z cm = 2.42 QSO J003 5-0918 was selected 
for our on-going survey for VMP DLAs (|Cooke et al.ll2010bl 1 
on the basis of its SDSS spectrum which shows a DLA at 
z a b s — 2.340 with no apparent associated metal lines. Such 
cases are highly suggestive of narrow and weak metal ab- 
sorption lines which are undetectable at the coarse resolu- 
tion and signal-to- noise ratio (S/N) of the discovery SDSS 
spectra. 

Follow-up observations of J0035— 0918 were made 
with the Magellan Echellette (MagE) spectrograph 
dMarshall et all 2008) on the Magellan n Clay telescope on 
the nights of 2008 December 30 and 31 in good conditions 
with 1 arcsec seeing. We took 2 x 2400 s exposures using 
the standard setup with lxl binning and a 1.0 arcsec slit 
giving a resolving power R = A/AA ~ 4100. Standard 
calibrations were taken following the recommendations by 
Simcoe, Hennawi & Williamo The data were reduced 
using a custom s et of IDL rout i nes w ritten by G. D. Becker 
and described in lBecker et al.l (|2006l ). The MagE spectrum 
confirmed the VMP DLA nature of the z abs = 2.3400972 
absorption system, by showing clear damping wings to 
the Lya absorption line (see top panel of Figure [1]) and 
unusually weak associated metal lines. 

Encouraged by these initial indications, we subse- 
quently observed J0035-0918 on the night of 2009 De- 
cember 9 with the High Resolution Echelle Spectrograph 



2 see http://web.mit.edu/rsimcoe/www/MagE/mage.html 



|Vogt et al.l 1 199^ ) on the Keck I telescope under good con- 
ditions with sub-arcsecond seeing, for a total integration 
time of 16200 s, divided into 6 exposures of 2700 s, result- 
ing in a signal-to-noise ratio near 4500 A of S/N ~ 18. 
We used the 1.148 arcsec wide slit which, when uniformly 
illuminated, provides a nominal spectral resolution R = 
A/AA = 37000. From our spectra, we measure R ~ 41000 
which corresponds to a velocity full-width at half maximum 
FWHM = 7.3 km s~ , sampled by ~ 3 pixels. We employed 
the UV cross-disperser which covers the wavelength range 
3100-6000 A with ~ 70 A-wide gaps near 4000 A and 5000 A. 

The HIRES spectra were reduced with the makee data 
reduction pipeline developed by Tom Barlow, which per- 
forms the usual steps of flat-fielding, order tracing, back- 
ground subtraction, and extraction of the final 1-D spec- 
trum. The data were wavelength calibrated by reference to 
the spectrum of a ThAr lamp, and mapped onto a vacuum 
heliocentric wavelength scale. The extracted spectra were 
merged and then normalized by dividing out the QSO con- 
tinuum and emission lines, using the software package uves 
POPLER, maintained by Michael Murphj[f|. Following this 
step, all available metal absorption lines associated with the 
DLA were extracted in ±150 km s _1 windows around the 
pixel with highest optical depth. A further fine correction to 
the continuum was then applied if necessary. 



3 PROFILE FITTING AND ABUNDANCE 
ANALYSIS 

Table [TJ lists all the metal absorption lines in the z a b s = 
2.3400972 DLA detected in our HIRES spectrum of 
J0035— 0918, together with a few interesting transitions 
which are below the detection limit of the data. For each line, 
we give the rest-frame wavelength and oscillator strength 
that we adopted for this work, as well as the measured equiv- 
alent width and associated random and systematic errors. 
The latter were determined by shifting the continuum by 
±a/^/n (where a is the 1<t error spectrum, and n is the 
number of independent resolution elements over which the 
equivalent width integration was carried out), and recalcu- 
lating the equivalent width. For the undetected transitions, 
we quote the 3a limiting rest-frame equivalent width, us- 
ing as a guide the profile of the weakest absorption line, 
Fell A1608, which we detect at the 5a confidence limit (see 
Table[TJ. Examples of metal absorption lines are reproduced 
in Fig. □ 

As can be appreciated from inspection of Table [1] and 
Fig-Hi the metal lines in this DLA are very weak and narrow, 
with equivalent widths Wo < 55 mA (the strongest line being 
CllA1334), and with the absorption taking place in a single 
velocity component with FWHM < 10 km s _1 . The weakness 
of the absorption limits our detection to the intrinsically 
most abundant elements of the periodic table, C, N, O, Al, 
Si, and Fe; on the other hand, the wide wavelength coverage 
of the echelle spectra, which reach well into the far-UV, gives 
access to several transitions of differing /-values for most of 
these elements. 



3 uves poplbr is available from 

http: / /astronomy.swin. edu.au/ ~mmurphy/UVES_popler 
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Figure 1. Selected absorption lines in the 2 a b s = 2.3400972 DLA towards J0035— 0918. The data are shown with black histograms, 
while the red continuous lines are model fits to the line profiles. Top panel: Portion of the MagE spectrum of J0035— 0918 encompassing 
the damped Lyo line, together with the theoretical Voigt profile (red line) for a neutral hydrogen column density log[A r (Hl)/cm -2 ] = 
20.55 zfc 0.1. The remaining panels display portions of the HIRES spectrum near metal lines of interest, together with model profiles 
generated with VPFIT as described in Section 13.11 This DLA consists of a single absorption component at z a b s — 2.3400972 with a small 
velocity dispersion, b = \f2a = 2.4km s — 1 . The weak absorption feature centred at +26 km s _1 in the Si II A1260.4221 panel is probably 
Fe II A1260.533 absorption in the DLA, although its strength is below our 3cr detection limit. The red wings of both C II A1036 and O I A988 
are blended with a weak Lycr forest line indicated by a continuous blue line. In all plots the y-axis scale is residual intensity. 
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Table 1. Metal lines in the 2 abs 
J0035-0918 



2.3400972 DLA towards 



Table 2. Ion column densities in the 
towards J0035-0918 



2.3400972 DLA 
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a Laboratory wavelengths and /-values from iMortonl l|2003h . 

b Rest-frame equivalent width. 

c Random error on the equivalent width Wo- 

d Systematic error on Wo from lcr uncertainty in the continuum 
placement. 

° 3cr upper limit on the rest-frame equivalent width. 



3.1 Column Densities 

We begin our abundance analysis by measuring the column 
density of neutral hydrogen. To this end, we used the MagE 
spectrum of the QSO, the relevant portion of which is repro- 
duced in the top panel of Fig.[T] because at z a b s = 2.3400972 
the damped Lya line falls on a gap between two of the CCDs 
in the HIRES detector mosaic. Even at the coarser resolu- 
tion of MagE (compared to HIRES), the broad damped Lya 
line is fully resolved and no loss of accuracy results in the 
derivation of iV(Hl). Fitting a Voigt profile to the line, we 
deduced log [iV(Hl)/cm~ 2 ] = 20.55±0.10; the corresponding 
theoretical Voigt profile is overlaid on the MagE spectrum 
in the top panel of Fig. [T] 

In the next step, we determined the Doppler parameter 
of the absorbing gas, b (km s _1 ), and the column density of 
the metal ions, iV(X) (cm -2 ), by fitting the corresponding 
line profiles with vpfitJj, which simultaneously fits multiple 
Voigt profiles to several atomic transitions, returning the 
values of N(X) and b, together with the associated errors, 
that minimize the y 2 between the data and the model. We 
tied the redshift and the Doppler parameter to be the same 
for all of the absorption lines listed in Table [TJ which is 
justified if the neutrals and first ions are kinematically asso- 
ciated with the same gas (we relax this assumption later). 
With these constraints, vpfit converged to a best-fitting 



Ion 


logAf(X)/cm- 2 


Hi 


20.55 ± 0.10 


Cn 


15.47 ± 0.15 


Ni 


13.51 ± 0.06 


Oi 


14.96 ± 0.08 


Aln 


11.73 ± 0.05 


Si n 


13.41 ± 0.04 


Sn 


< 13.08 


Fen 


12.98 ± 0.07 



model consisting of a single absorption component with red- 
shift z = 2.3400972 ± 0.0000008 and Doppler parameter 
b = 2.36 ± 0.08 km s . The corresponding column den- 
sities are listed in Table [2] Examples of the theoretical line 
profiles generated by vpfit are shown superimposed on the 
data in Fig. Q] 

The weakest feature in our data is SllA1259. When 
this line is included in the vpfit fitting procedure (see 
bottom right panel of Fig. [!}, we derive a column density 
logiV(Sn)/cm~ 2 = 13.08 ± 0.10. However, since this ab- 
sorption line is only significant at the ~ 4.5a level, we con- 
servatively consider the above value to be an upper limit to 
the column density of S II. 

With the usual assumption that the ions observed are 
the dominant stage of the corresponding elements in H I gas, 
so that corrections for unseen ion stages and/or the presence 
of ionized gas are negligible (we review this assumption in 
Section [3.21 below) . it is straightforward to deduce the abun- 
dances of the elements concerned by simply dividing the 
column densities in Table[2]b y iV(Hl). Compari son with the 
solar abundance scale of Asplund et al.l (|2009l ) then leads 
to the abundance pattern listed in Table and illustrated 
graphically in Fig. [3] 



3.2 Ionization Corrections 

In DLAs, it is usually assumed that the metals within the ab- 
sorbing Hi gas reside in a single dominant ionization stage, 
Xn, so that the total abundance of an element is given by 



[X/H] = [Xn/Hi] +IC(X) 



(1) 



vpfit is available from http: / /www. ast.cam.ac.uk/ ~rfc/ vpfit. html 



where the ionization correction, IC(X), is typically negligi- 
ble. In general, it is safe to assume IC(X) ~ 0.0 for gas with 
high iV(H l), since the gas is s elf-shielded from ionizing radi- 
ation (e.g. lVladilo et all lioOl). Of course, if the gas does not 
reside in a single dominant ionization stage, or some amount 
of X N is associated with H n gas, we may respectively under- 
or over-estimate the abundance of a given element [positive 
or negative IC(X)]. 

To gauge the extent of such c orrections, we used the 
CLOUDY photoionization software (jFerland et all Il99&1 ) to 
model the DLA as a slab of gas with uniform density 
in the range —3 <_log[n(H)/cm -3 ] < 3, exposed to the 
lHaardt fc Madaul |200ll) metagalactic ionizing background 
and the cosmic microwave background, both at the red- 
shift of the DLA. A dopting the solar abundance scale of 
lAsplund et all l|2009tl . we globally scaled the metals to log 
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Table 3. Element abundances in the 
towards J0035-0918 



2.3400972 DLA 



Element 


log 


c(X) DLA lo 


gc(X)^ b 


[X/HJdla 




-1 





C 




6.92 


8.43 


-1.51 




-1 


5 


N 




4.96 


7.83 


-2.87 













6.41 


8.69 


-2.28 




-2 





Al 




3.18 


6.44 


-3.26 


\ 






Si 




4.86 


7.51 


-2.65 


X 


-2 


5 


s 




< 4.53 


7.14 


^ -2.61 








Fe 




4.43 


7.47 


-3.04 





















-3 


a loge(X) = 12- 


Ho 


rN(X)/N(R). 








-3 


5 



b Asplund et al. (2009). 
c [X/H] DLA = lo ge (X) DLA 
Table 2. 



log6(X) Q , 



with errors as listed in 



Zdla I Zq = —2.75 (an approximate average metallicity — 
see Fig. [2}. Once the column density of the DLA was 
reached, we stopped the simulations and output the result- 
ing ion column densities. Using these computed column den- 
sities, we are then able to derive the ionization corrections 
from Eq. Q] 

The results of our CLOUDY simulations are shown in 
Fig. [3] The ionization corrections appropriate to the DLA 
under investigation depend on the volume density of the gas, 
which can be inferred from the ratio of successive ion stages 
(see right panel of Fig. [3J. The only element for which we 
have this information is nitrogen and, even then, we can 
only derive a 3cr upper limit to the N II column density from 
the 3a upper limit on the rest frame equivalent width of the 
undetected N II A1084 line (see Table[l|, leading to the upper 
limit log jV(Nn)/iV(Ni) < -0.91. This affords a lower limit 
on the volume density of log[n(H)/cm -3 ] ^ —1.0 (see right- 
hand panel of Fig. [3}. Referring now to the left-hand panel 
of Fig. [31 it can then be seen that, when the gas density 
is greater than 0.1cm -3 , the ionization corrections for the 
ions of interest here are indeed small, IC ,$ 0.1 dex. These 
values are comparable to the uncertainties in the ion column 
densities (see Table [2} which justifies our assumption that 
[X/H] ~ [Xn/Hi], and that ionization corrections are not a 
serious concern. 

It is worth noting, in passing, that even in the ab- 
sence of several density-sensitive ion ratios, low values of 
n(H) are in general unlikely in these metal-poor DLAs 
with such simple velocity structure. The column density 
log[AT(Hl)/cm~ 2 ] = 20.55 implies a linear size / > 1 kpc 
along the line of sight if log[n(H)/cm -3 ] < —1. It seems 
unlikely that the structures giving rise to DLAs with prop- 
erties similar to that considered here could have such large 
physical extent while maintaining very quiescent kinematics, 
with internal velocity dispersions of only a few km s -1 . 



4 A MASQUERADING CARBON 
ENHANCEMENT? 

Returning to Fig. H we note that with [O/H] = -2.28 
and [Fe/H]= -3.04 the z = 2.3400972 DLA in line 
to J0035— 09 1 8 is among the mo s t me tal-poor known 
|Pettini et al.l 120081 ; IPenprase et all l2010l l. However, the 



Al 



Si 



Fe 



Figure 2. Element abundances in the z abs = 2.3400972 DLA 
to wards J0035— 0918 co mpared to the solar values as compiled 
by lAsplund et all J2009T ) . The height of each box represents the 
uncertainty in each element abundance; the 3<j upper limit for S 
is indicated by the bar and arrow. 



most striking feature of the abundance pattern in Fig. [2] is 
the overabundance of carbon relative to all other elements. 
Thus, for example, [C/O] = +0.77, which implies an over- 
abundance of carbon relative to oxygen by a factor of ~ 6; 
similarly, [C/Fe] = +1.53, or ~ 35 times solar! Given such 
extreme values, it is important to consider what factors, if 
any, may have led us to spuriously overestimate the abun- 
dance of carbon. To this end, we performed several tests, 
which we now discuss. The upshot of these tests is that we 
could not find any plausible reason why the abundance of 
carbon would have been significantly overestimated. Thus, 
readers who are less interested in the details of these checks 
may skip the rest of this section and move straight to the 
discussion in Section [5] 



4.1 Incorrect background subtraction? 

First we considered the possibility that the background may 
have been oversubtracted in the proximity of the C 11 absorp- 
tion lines; if the zero level had been incorrectly determined, 
this may lead us to overestimate the apparent optical depth 
of the absorption. While we cannot categorically rule out 
this possibility, as we do not have any independent mea- 
sures of the zero level (such as saturated absorption lines 
with flat cores) in the immediate vicinity of the Cll lines, 
we note the following: (i) there are a number of strongly sat- 
urated Lya absorption lines in the Lya forest (shortwards of 
Aobs = 4158 A) and none of them show a systematic offset 
of their cores from zero residual intensity; (ii) the appar- 
ent optical depths of the two C n transitions covered by our 
HIRES spectrum, A1334 and A1036, are mutually consistent 
(see Fig. [TJ; thus, if the background level were incorrect, it 
would have to have been oversubtracted by the same frac- 
tional amount of the QSO continuum at these two wave- 
lengths, separated by 995 A in the observed spectrum (at 
z abs = 2.3400972). We also inspected the raw 2-D HIRES 
frames in the region of the two best observed C 11 and O I 
lines, A1334 and A1302 and found that they both fall close 
to the peak of the echelle blaze function, near the centre of 
the HIRES detector, where the data are of the highest S/N. 
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Figure 3. Left panel: Ionization corrections, as denned by Eq. [T] for elements of interest (symbols connected by a solid line), as a 
function of the volume density of the gas in the z^ s = 2.3400972 DLA. Right panel: Ratio of the computed column densities for the 
successive ion stages of nitrogen (symbols connected by a solid line), plotted as a function of the volume density. The solid horizontal 
line indicates the upper limit we measure for log N(N Il)/jV(N i), whilst the vertical line with an arrow indicates the corresponding lower 
limit on the volume density of the gas, log[n(H)/cm~ 3 ] ^ —1.0. Both plots were produced with CLOUDY photoionization calculations; see 
the text for further details. 



We conclude that the recorded excess of C is unlikely to be 
an artifact of the data reduction process. 

4.2 Profile fitting 1: gas kinematics 

Next we looked critically at the profile fitting procedure. 
Our vpfit modelling described in Section \3 . 1 1 assumed that 
the neutrals and singly ionized species arise from gas with 
the same Doppler parameter. We relaxed this constraint by 
fitting separately the absorption lines of the first ions (two 
C II, five Si n lines and one line each of Al n and Fe n) and the 
neutrals (five Ol and two Nl lines). In this case, vpfit con- 
verged to best fitting Doppler parameters of b — 2.30 ± 0.09 
for the first ions and b = 2.6 ± 0.2 for the neutrals. With 
these parameters the C II column density is higher by 0.1 dex 
than the value listed in Table [2] (that is, the value obtained 
by tying b to be the same for all species), while N(Oi) is 
lower by 0.1 dex. The column densities of other elements re- 
main essentially unchanged. In other words, kinematically 
decoupling neutrals and first ions has the net effect of fur- 
ther increasing the [C/O] overabundance relative to the val- 
ues in Table [3] and Fig. [2] Note also that the above values 
of b are still consistent within the errors with the best- 
fitting b = 2.36 ± 0.08 deduced in Section 13.11 by tying 
all the absorption lines together. Finally, fitting only the 
C II lines, without reference to any other transition, yielded 
log[iV(Cll)/cm" 2 ] = 15.42 ± 0.59, which is only 0.05 dex 
lower than the value listed in Table [2l albeit with a larger 
error (but with much the same redshift and Doppler param- 
eter: z abs = 2.340096 ± 0.000002 and b = 2.4 ± 0.3km s" 1 ). 

4.3 Profile fitting 2: instrumental resolution 

In order to compare theoretical and observed absorp- 
tion line profiles, vpfit requires knowledge of the in- 
strumental broadening function, normally assumed to be 
a Gaussian with FWHM corresponding to the spectral 
resolution. The nominal resolution of spectra recorded 
through the 1.148 arcsec wide entrance slit of HIRES em- 
ployed in our observations is FWHM = 8.1km s _1 (see 



http://www2.keck.hawaii.edu/inst/hires/), but this applies 
to a uniformly illuminated slit. Since the seeing was con- 
sistently better than 1.148 arcsec during our observations of 
J0035— 0918, it is likely that the actual resolution of our data 
is somewhat better than the nominal value. The narrowest 
features in the spectrum of J0035— 0918 are the metal ab- 
sorption lines of the z a t, s = 2.3400972 DLA themselves. We 
therefore estimated the true spectral resolution — measured 
by the Doppler parameter &i ns tr = 0.6006 FWHM — by vary- 
ing bmstr in small steps from 5.0 to 3.7 km s _1 (i.e. between 
FWHM = 8.3 and 6.2 km s" 1 ), fitting all of the absorption 
lines as described in Section \S. II and minimizing the x 2 De " 
tween the model and observed profiles as a function of fcinstr- 
This procedure gave a well-defined bi ns tr,min = 4.4 km s , 
corresponding to FWHM = 7.3 km s _1 , or R = 41 000, which 
is the value used in all the profile fitting described in Sec- 
tions O and |4l 

It is obviously important to test the sensitivity of our 
results to the value of &j nstr adopted. Within the range of 
values tested, &i ns tr = 5.0-3.7 km s _1 , the column density 
of C II and the corresponding element ratios changed by no 
more than ±0.1 dex compared to the values in Table[5]which 
refer to the case foinstr.mm = 4.4 km s _1 . On the basis of the 
tests carried out in this and the preceding subsection, we 
conclude that the overabundance of carbon is not an artifact 
of the profile fitting analysis of the absorption lines. 

4.4 Monte Carlo simulations 

The z = 2.3400972 DLA in J0035-0918 is unique so far 
in showing such a marked overabundance of carbon. More 
typically, DLAs with [O/H] < -2 have [C/O] ~ -0.2 and 
at most [C/Q1 ~ 0.0 [i.e . (C/O) = (C/Q) B ; IPettini et al.l 
l200St ICooke et alJ|2uTobl |. It is therefore worthwhile consid- 
ering to what extent a solar C/O ratio is excluded by our 
data. We illustrate this test in Fig. [4j where we compare the 
line profiles for the two C II absorption lines generated with: 
(upper panels) the best-fitting model returned by vpfit, and 
(lower panels) a model in which N(C n) has artificially been 
fixed at the column density corresponding to [C/O] =0.0 



© 2010 RAS, MNRAS 000.[TJfT3l 



8 Cooke et al. 



(a) [C/0] = 0.77 
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Figure 4. Comparison between observed C II line profiles (black histograms) and theoretical profiles computed with VPFIT (red continuous 
lines). Below each plot we also show the difference between computed and observed residual intensities (red lines) compared with the 
la error spectrum of the data (grey area), (a): best fitting model as described in Section |3. 1 1 with the log[W(C ll)/cm — 2 ] = f5.47 as in 
Table [2] (b): best fitting model obtained by forcing log[JV(C ll)/cm~ 2 ] = 14.70 so that [C/O] = 0.0 (see text for further details). This 
model provides a poorer fit to both Cn absorption lines, as demonstrated by the mismatch between the model subtracted from the data 
and the la error spectrum (bottom panels). 



(that is log[iV(Cn)/cm~ 2 ] = 14.70 as opposed to 15.47 as 
in Table[2j). While the former results in a xVdof = 13.9/22, 
the latter fits are worse, with x 2 /dof = 86.6/22. The higher 
value of x 2 is reflected by the higher residuals, both in the 
core and the wings of the C II A1334 line in particular, as can 
be appreciated from close examination of Fig. [4] The nois- 
ier CnA1036 line is less instructive in this context, partly 
because of blending with Lya forest lines. 

It may be surprising to some to see what subtle changes 
in the line profiles result from changing N(C n) by a factor 
of ~ 6. The reason is that the two C n absorption lines are 
close to saturation and, as they approach the flat part of 
the curve of growth, their equivalent widths only increase 
slowly with rising column density. In the present case, the 
curve of growth for the metal lines in the z a b s = 2.3400972 
DLA towards J0035-0918 (reproduced in Fig. © is well- 
constrained by the relative strengths of six O I and five Si II 
transitions of widely differing /-values (see Table [T]). Since 
the column densities of both O I and Fe n are fixed by op- 
tically thin transitions (i.e. independent of the Doppler pa- 
rameter), the ratios [C/O] and [C/Fe] depend only on the 



C n column density. In the right panel of Fig. [5] we com- 
pare the measured equivalent width of CllA1334 (indicated 
by the blue hatched region) with the values expected if 
[C/O] =0.0 (red), [C/Fe] =+1.0 (blue), and [C/Fe] = +1.53 
(green). Only in this last case ([C/Fe] = +1.53) do we re- 
cover the measured Wo(C II A1334), whereas the strength of 
the line is underpredicted in the other two cases illustrated. 
The partially blended CllA1036 line is not instructive in 
this equivalent width test which, by its nature, is less sen- 
sitive than the pixel-by-pixel line profile analysis illustrated 
in Fig. II 

Figures [4] and [5] offer a clearer view of the carbon- 
enhanced nature of this DLA, while at the same time high- 
lighting the subtle differences in the absorption line profiles 
when lines are on the flat part of the curve of growth. Thus, 
it is reasonable to question whether the differences between 
the line profiles computed with [C/O] = 0.0 and the observed 
profiles, most evident for C II A1334 in the bottom right-hand 
set of panels in Fig. 21 may be simply due to statistical fluc- 
tuations. We tested this possibility with a Monte Carlo-type 
approach. We synthesized 100 pairs of CllA1036, A1334 line 
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log [ N/cm- 2 ] log [ N(C II) /cm- 2 ] 

Figure 5. Left panel: Curve of growth for the z a b s — 

2.3400972 DLA in line to J0035-0918, with Doppler parameter b = 2.36 ± 
0.08 km s — 1 . The dashed box shows the region reproduced on an expanded scale in the right panel. Right panel: Equivalent width 
comparison. The continuous black line is a portion of the curve of growth for Doppler parameter b = 2.36km s — x , with the short-dash 
lines indicating the uncertainty Sb = ±0.08 km s _1 . The red, blue, and green horizontal bands show the corresponding ranges in equivalent 
width Wo(CllA1334) appropriate to element ratios [C/O]=0.0, [C/Fe] = +1.0, and [C/Fe] =4-1.53 respectively. Only the last of these 
(green band) matches the measured equivalent width of CllA1334, indicated by the blue hatched horizontal band. A strict lower limit 
[C/Fe] ^ +1.25 (vertical long-dash line with right-pointing arrow) is obtained by considering the lower bound on Wo(CllA1334) and 
upper bound on b. 



profiles, all corresponding to [C/O] = 0.0 and with the best 
fitting Doppler parameter (b = 2.36 ± 0.08). We perturbed 
these line profiles with Gaussian distributed errors from the 
observed la error spectrum, and then used vpfit to fit each 
of these 100 pairs of Cll lines (together with the observed 
line profiles of the other ions) and deduce 100 new values of 
N(Cll). Using the simultaneously derived Ol column den- 
sity (which is essentially unchanged since it is constrained by 
weaker transitions), we calculated the [C/O] ratio for all of 
our realizations to determine how often one would expect to 
observe [C/O] ^ 0.77 if the true [C/O] ratio were solar. The 
results of this Monte Carlo approach are shown in Fig. [5] 
In none of the simulations was a value of [C/O] as high as 
+0.77 recovered. The distribution of recovered values is sym- 
metric about the input value [C/O] =0.0 with a dispersion 
a = 0.087. Thus, if the [C/O] ratio in the z = 2.3400972 
DLA were indeed solar (which in itself is higher or as high 
as observed in any other DLA up to now), our data would 
constitute a ~ 9<r deviation, which is very unlikely indeed. 

In conclusion, none of the tests we have performed can 
explain away the C overabundance we have uncovered. The 
most straightforward interpretation of the data is that we 
have identified a damped Lya system with a marked over- 
abundance of C, comparable to that so far found only in 
carbon-enhanced metal-poor stars of the Milky Way. 



5 A CARBON-ENHANCED METAL-POOR 
DLA 

As already mentioned, this is the first example of a DLA 
meeting the conventional definition of CEMP stars, [C/Fe] > 
+ 1.0. While it is t he only such case in the survey by 
ICooke et all l|2010bD . which includes a dozen DLAs with 
[Fe/H] < —2, it is nevertheless very hard to draw conclusions 
regarding the frequency of CEMP DLAs. The reason stems 
from the difficulties associated with measuring the Cll col- 
umn density; the CllA1036 and A1334 absorption lines are 
almost always saturated in damped Lya absorption systems, 
even at low metallicities. CEMP DLAs may thus be easily 
overlooked. 




-0.3 -0.2 -0.1 0.1 0.2 0.3 



[C/O] 

Figure 6. Histogram of the values of [C/O] recovered from 100 
Monte Carlo realizations with input [C/O] = 0.0 in which the syn- 
thetic C II A1334, A1036 line profiles were perturbed with a ran- 
dom realization of the error spectrum. The overlaid Gaussian fit 
to the results has a dispersion a = 0.087, which implies that our 
measured value [C/O] = +0.77 would be a ~ 9<r random fluctua- 
tion if the [C/O] ratio were solar in the z = 2.3400972 DLA. 



The discovery of a gaseous, high redshift, counterpart 
to CEMP stars of the Galactic halo is a breakthrough with 
potentially very important consequences. As mentioned in 
the Introduction, interpretations still differ as to the ori- 
gin of the carbon enhancement in very metal-poor stars, 
but the favoured explanation involves mass transfer from 
an AGB companion, at least for the members of the CEMP 
class which also exhibit enhancements in s-process elements. 
The situation is less clear-cut for the CEMP stars with 
no such enhancement, the CEMP-no stars. On the one 
hand, there may be a c ontinuous link betwe en this class 
and the CEMP-s class (|Masseron et al.ll20ldh . so that in 
some CEMP-no stars the carbon enhancement may also be 
due to mass transfer from an unseen companion. On the 
other hand, the binany transfer scenario cannot readily ex- 
plain why the fraction of all metal-poor stars that exhibit 
a carbon-enhancement in creases with decreasing metallicity 
IIBeers fc Christliebll2005h . In fact, the only physical model 
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so far put forward that can explain this increased fraction re- 
quires a high carbon yield from metal-free stars to efficiently 
cool low metallicity gas a nd drive the transition from Pop II I 
to Pop II star formation (|Frebel. Johnson. fe Brommll2007l ). 

While CEMP-no stars in the Galactic halo may have 
more than one origin, the same ambiguities do not apply to 
the carbon enhancement in the DLA presented here, where 
we have an entire cloud overabundant in carbon (below we 
place limits to the total mass of carbon involved). Rather, 
what we presumably see in the DLA is the initial chemi- 
cal composition of the gas from which CEMP stars would 
subsequently form. In this case, the pattern of element abun- 
dances revealed gives strong clues to the nature of the earlier 
generation of stars (and supernovae) which seeded the DLA 
cloud with its metals. 

We explore these ideas further by comparing the abun- 
dance pattern of the DLA with those of halo stars of simi- 
lar metallicity and with calculations of element yields from 
metal-free and metal-poor stars. In this endeavour, we are 
limited by the small number of elements whose abundances 
we have been able to measure, but this is an inevitable con- 
sequence of the low metallicity of this DLA. While it would 
be instructive to know the abundances of other Fe-peak ele- 
ments, such as Cr and Zn for example, the absorption lines 
of these species would require S/N > 500 for 5a detections^ 
Realistically, only Mg and S are within reach of dedicated 
future observations, and there is certainly no hope of mea- 
suring directly the abundances of neutron capture elements 
in the DLA. 



5.1 Comparison with Stellar Abundances 

In order to compare the relative abundances of the ele- 
ments measured in the DLA with the corresponding values 
in Galactic halo stars, we median-combined the abundances 
of all CEMP stars with Fe/H within a factor of two of the 
D LA (-3 . 34 < [Fe/H] < -2.74) from the recent compilation 
of iFrebell l|2010h . We then calculated the dispersion in each 
abundance using the Interactive Data Language routine RO- 
BUST.SlGM^fl which determines the median absolute devia- 
tion (unaffected by outliers) of a set of measurements, and 
then appropriately weights the data to provide a robust esti- 
mate of the sample dispersio n (|Hoaglin. Mosteller. fe Tukevl 
1 19831 . cf. lCooke et~aDl2010ah . 

Figure illustrates three comparisons. In the top panel, 
the DLA abundance pattern (solid black boxes) is shown 
together with that of the 'average' population of all CEMP 
stars (open magenta boxes), while in the middle panel the 
comparison is restricted to CEMP-no stars. Clearly CEMP 
stars, of both flavours, are a heterogeneous population, ex- 
hibiting wide ranges in the ratios of the elements considered 
here. On the basis of the limited evidence available, it ap- 
pears that the C/Fe ratio in our DLA is a better match 
to the values seen in the subset of CEMP-no stars than to 
the more extreme values encountered in some CEMP stars. 
Other elements ratios are less straightforward to interpret. 

Nitrogen, for example, is evidently less abundant in the 
DLA compared to CEMP stars by ~ 2 orders of magnitude! 



5 Assuming [Cr/ Fe] = [Zn/Fe] = 0.0 . 

B Available from http://idlastro.gsfc.nasa.gov/homepage.html 
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Figure 7. Comparison of element abundances in the z^ s = 
2.3400972 DLA (filled black boxes) and in Galactic halo stars 
with an Fe abundance within a factor of two of the DLA (open 
magenta boxes). The numbers below the element labels indicate 
the number of stars that contributed to the determination of the 
'typical' stellar abundances, and the heights of the magenta boxes 
reflect the dispersion of each set of measurements. Top Panel: 
comparison with all CEMP stars that have -3.34 sg [Fe/H] < 
—2.74. Middle Panel: comparison with CEMP-no stars that have 
—3.34 ^ [Fe/H] —2.74. For this case, the oxygen abundance 
of a single CEMP-no star is shown by the open circle. Lower 
Panel: comparing the DLA abundance pattern with the stellar 
abundance patterns of HE 0143-0441 (red symbol s; a CEMP- 
s sta r with [Fe/H] = -2.21, [Ba/Fe] = +0.62 from ICohen et al.l 
I2004F) and BD+44°493 (green symb ols; a CEMP-n o star with 
[Fe/H] = -3.73, [Ba/Fe] = -0.55 from llto et aDbOOSh . Note that 
in this last panel, we have plotted [X/Fe] as opposed to [X/H]. In 
all panels the dashed line represents the solar abundance. 



However, it is difficult to draw firm conclusions from this dis- 
crepancy, because the abundance of N is notoriously difficult 
to measure in metal-poor stars where it relies on the analysis 
of molecular band s of NH or CN , each with its own disad- 
vantages (see e.g. lAsplundll2005t ). The largest uncertainty 
arises from the assumed ID (rather than 3D) model atmo- 
spheres; correcting for 3D effects can change then deduced 
values of N/H by ~ —0.5 dex (Note that similar corrections 
also apply to C/H when derived from molecular bands of CH 
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or C^ ICollet. Asplund. fc Tramped ach 2007;). Furthermore, 
the photospheric N/H we measure today is unlikely to re- 
flect the N abundance of the gas cloud that gave birth to the 
star; during the star's evolution it will become self-polluted 
by its own nucleosynthesis (e.g. through CN cycling or rota- 
tional mixing), which is likely emphasised by the high seed 
C abundance. Thus, a comparison between stellar and DLA 
nitrogen abundances is not very instructive. 

Comparing the oxygen abundances in CEMP stars and 
the DLA is similarly problematic. The statistics are very 
poor (see Fig. [7|) and it is well known that different spectral 
features used to ded uce O/H in the metal-poo r regime give 
discordant answers jGarcfa Perez et al l 120061 ). Significant 
departures from local thermodynamic equilibrium (LTE) be- 
come important when measuring the O abundance from 
the infrared triplet at 777 nm in the metal-poor regime 
([Fe/H] < —2.5), leading t o non-LTE correction s of the or- 
der of -0.5 to -1.0 dex (|Fabbian et al.l 1200901 ). Similarly, 
the oxygen abundance determined from the UV OH lines 
is subject to l arge 3D correct ions of up to —0.9 dex when 
[Fe/H] sC -3.0 (|Asplundll20"05l '). Clearly, it will be important 
to reexamine this issue once samples of CEMP stars have 
grown to include more cases where [O/Fe] has been mea- 
sured with confidence. 

Finally, Fig. [7] demonstrates that there is a reasonably 
good match between CEMP stars and the DLA in the el- 
ements Al and Si, particularly for the CEMP-no stars. In 
this context, we point out that significant (positive) non- 
LTE corrections to A l/H may apply for stars in this metal- 
licity regime (see e.g. lAsplundl2005h . To address this issue, 
we first identified those CEMP stars in the iFrebell (|20ld ) 
compilation whose Al abundances were derived assuming 
LT E. We then applied non- LTE corrections to Al based on 
the lAndrievskv et alj |2008t ) estimates, which amounted to 
a typical correction of the order of +0.6. The Al abundances 
for the 'typical' CEMP and CEMP-no star presented in 
Fig. [7] should thus correspond to the (approximate) non- 
LTE values. 

Finally, in the lower panel of Fig. [7] we compare the 
element ratios (relative to Fe) in the DLA with the values 
me asured for two CEMP stars selected from the compilation 
bv lFrebell (|2010T ) because they most closely match the rel- 
ative abundances in the DLA. It is certainly plausible that 
stars forming out of this DLA gas would share many chem- 
ical similarities - at least for the elements considered here - 
with these two stars of the Milky Way Galactic halo. 



5.2 Comparison with Model Yields for Metal- free 
Stars 

With a metallicity [Fe/H] of only 1/1000 solar, it is con- 
ceivable that the chemical composition of the DLA gas re- 
flects the element yields of only a few prior generations of 
stars. It is thus of interest to compare the abundance pat- 
tern in Fig. [2] with calculations of nucleosynthesis by metal- 
free, or low metallicity, stars. Calculations of yields from 
metal-free stars (commonly referred to as Pop III stars), 
aimed in particular at interpreting the element ratios seen 
in very metal-poor stars, have foc used on models wh i ch in- 
clude 'mixing and fallback' (e.g. lUmeda fc Nomotol 12001 
|2003I :Ih" eger fc Wooslevllioioh . These scenarios involve core- 
collapse supernovae where the elements synthesized in the 



inner regions of the star are mixed by some process (several 
possibilities have been put forward); a fraction of the mixed 
material subsequently falls back onto the central remnant 
while the rest is ejected into interstellar space. 

Such models are illustrative, rather than predictive, 
since the parameters describing the boundaries of the mixing 
region and the fraction of the mixed material which is ejected 
cannot be derived from first principles, and are instead pa- 
rameterised and suitably adjusted to fit the observed stellar 
abundances (e.g. iTominaga et al.l [2007^ . A string of recent 
models that em ploy a more physically motivated prescrip- 
tion of fallback l|Heger fc Wooslevll2010f ). have also investi- 
gated the effects of mixing due to the Ray leigh- Taylor insta- 
bility (|Joggerst. Wooslev. fc H cgcr 200 9]) together wit h ro- 
tati onally induced mixing in 2D djoggerst et alJfeoiOal ) and 
3D (jjoggerst, Almgren. fc Wooslevl l2010bl ). By and large, 
a general feature of these models is that an abundance 
pattern similar to that found here in the DLA, including 
the carbon enhancement and the marked odd-even effect 
(see Fig. [2|, can be reproduced with low energy explosions, 
su ch as those giving rise to the faint SN branch (see Fig. 1 
of ITominaga et al.ll2007l . and Fig. 11 of iHeger fc Wooslevl 
120101 ) , and a moderate degree of mixing and fallba ck. For ex- 
ample , in the mixing-fallback model considered bv lTominagal 
(2010). enrichment with the elements ejected by a single 
25 Mq star of zero metallicity, exploding as a core-collapse 
supernova with an energy _Esn ~ 10 51 ergs, provides a re- 
markably good fit to the relative element abundances mea- 
sured in the DLA (Kobayashi et al. in prep). 

It remains to be established whether a zero initial metal- 
licity (i.e. Pop III) is actually required to match the ob- 
servations, or whether a similar abundance pattern could 
also be reproduced by models with low, but non-zero, ini- 
tial metallicity. An alternative to the mixing and fallback 
model has been proposed by the Geneva group whose work 
places more emphasis on the effects of rotation to provide 
the mixin g and the trigger for mass loss in very metal-poor 
stars (e.g.TMevnet, Ekstrom, fc Maederll2006l ; lHirschi|[2007l ; 
iMevnet et al.|[2Q10l ) . Whilst these models are physically well- 
motivated, the rather high nitrogen yield predicted by some 
models contrasts with the relatively low [N/C] and [N/O] 
observed in this DLA. This may in turn place interesting 
contraints on the rotation velocities of the metal-free stars 
that may have seeded the DLA with its metals. Only re- 
cently has rotation been included in zero- metallicity pro- 
genito rs under the mixing-fallback scenario jjoggerst et al.l 
l2010al ). It should be possible to assess better the effects of 
rotation on the yield from zero-metallicity progenitors as 
these stellar models improve and as the samples of CEMP 
DLAs grow. 



5.3 How Many Supernovae? 

It is truly intriguing that an entire cloud of neutral gas is 
so highly enhanced in carbon and shows such a clear-cut 
odd-even effect, if both are indeed characteristic of Pop III 
supernova yields. In this context, it is of interest to estimate 
the mass of carbon involved. With the assumption of spheri- 
cal symmetry, the mass of singly ionized carbon in the DLA 
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is given by: 
M(Cll) = 12M(Hl) 

M(Cll) = 2 7rm H 



iV(Cll) 
AT (Hi) 
iV(Hi) 2 7V(Cii) 



n(H) 2 



M(Cn) sC 200 



"(H) 
0.1cm" 



(2) 
(3) 
(4) 



where run is the mass of a hydrogen atom and n(H) is the 
volume density (cm -3 ). The corresponding mass of neutral 
gas in the DLA is: 



M lo ^ Af(Hl) 3 



(5) 
(6) 



Since C is mostly singly-ionized in DLAs, eq. (4) im- 
plies an upper limit of 200 Mq to the total mass of 12 C 
in the DLA. However, the inverse square dependence of 
M(C n) on the gas density allows substantially lower val- 
ues of M(Cll). In Section [3.21 it was shown that our lower 
limit n(H) ^ 0.1 cm~ 3 was based on the non-detection of 
Nil absorption, and that higher densities are likely given 
the very low velocity dispersion of the gas. If, for example, 
n(H) ~ 1cm -3 , which would imply a linear extent of the 
DLA of 100 pc along the line of sight, then the implied mass 
of 12 C would be reduced to only 2 Mq. For comparison, the 
total mass of 12 C ejected b y the s ingle 25 Mq Pop III star 
in the model by iTominaeal J2010l ') is ~ 0.2M Q . The mod- 
els bv lHeeer fe Wooslevl < |201cF " anticipate a similar ejected 
mass of 12 C (~ 0.3 Mq) for a comparable explosion energy 
and progenitor mass. Thus, we may indeed be seeing the ele- 
ments synthesized by only a few supernovae in the chemical 
enrichment of the DLA considered here. 

In conclusion, we speculate that the z a bs = 2.3400972 
DLA in front of the QSO J0035-0918 may well be the much 
sought 'missing link' between the first, zero-metallicity, stars 
and the most metal-poor stars in the halo of our Galaxy. Its 
low metallicity of 1/1000 solar in Fe, coupled with an over- 
abundance of C and a marked odd-even effect in the relative 
abundances of the few elements that could be measured, are 
all consistent with the yields produced by models of Pop III 
stars which explode as core-collapse supernovae of relatively 
low energy. The mass of newly synthesized elements may be 
that produced by only a few such supernovae, depending on 
the unknown volume density of the gas in the DLA. In this 
scenario, the gas we see as a damped Lya system at high 
redshift may be the material from which a subsequent gen- 
eration of stars formed, with a chemical composition similar 
to that seen in CEMP-no stars of the Galactic halo. 



6 SUMMARY AND CONCLUSIONS 

In the course of our program to study the most metal- 
poor DLAs and, in particular, to measure the abundances 
of the CNO group of elements, we have uncovered a DLA, 
at z = 2.3400972 towards the z cm = 2.42 SDSS QSO 
J0035— 0918, whose chemical composition is consistent with 
that produced from exploding Population III stars. From the 



analysis of medium and high resolution spectra of this DLA 
we draw the following conclusions. 



(i) The metal absorption lines associated with the DLA are 
formed in quiescent gas, consisting of a single absorption 
component with a small Doppler parameter b — 2.4 km s . 

(ii) The metallicity of the DLA, as measured from Fe, is 
very low: [Fe/H] = -3.04 ±0.17, or ~ 1/1000 solar. 

(iii) The DLA exhibits a strong enhancement of the abun- 
dance of carbon relative to all other metals covered by our 
data: N, O, Al, Si, S and Fe. We measure [C/Fe] = +1.53, a 
factor of ~ 20 greater than observed in any other DLA up 
to now. Adopting the defining criterion for carbon-enhanced 
metal-poor stars in the Galactic halo, [C/Fe] > +1.0, this 
is the first example of an analogous carbon-enhanced DLA. 
We also deduce [C/O] = +0.77, whereas in all other DLAs 
with [Fe/H] < -2 studied up to now [C/O] < 0.0. 

(iv) The DLA also exhibits a clear odd-even effect, which 
implies a low neutron excess, and hence presumably low 
abundances of neutron-capture elements. When its chemi- 
cal composition is compared with that of Galactic carbon- 
enhanced metal-poor stars which do not exhibit an excess of 
neutron-capture elements, a good match is found for some 
element ratios (C/Fe, Al/Fe and Si/Fe). N and O are signif- 
icantly less abundant, compared to Fe, in the DLA than in 
most CEMP stars, but it is difficult to draw definite conclu- 
sions for these two elements which are notoriously difficult 
to measure in very metal-poor stars. 

(v) The abundance pattern we observe for this DLA is con- 
sistent with enrichment from a population of ~ 20-50 Mq 
metal-free, or extremely metal-poor, stars that ended their 
lives as core-collapse supernovae with modest explosion en- 
ergies. We estimate the total mass of 12 C within the DLA to 
be 200 Mq . This upper limit could be constrained further 
with higher signal-to-noise ratio spectra that would permit 
a measure of the column density of N n and other successive 
ion stages whose ratios are density dependent. The steep de- 
pendence of M( 12 C) on the gas density allows the possibility 
that we may be seeing the chemical enrichment produced by 
only a few prior supernovae. 

(vi) We speculate that the gas in this DLA may be the 
'missing link' between the yields of Population III stars, and 
their later incorporation in the CEMP-no class of carbon- 
enhanced metal-poor stars. We note, however, that the 
carbon-enhancement in some CEMP-no stars could also be 
produced by other means. Long term radial velocity moni- 
toring of CEMP-no stars will confirm or deny their associa- 
tion with a now extinct companion. 

The results presented here emphasise the importance 
of further observations of DLAs with [Fe/H] ^ —3 to com- 
plement the work being carried out on the most metal-poor 
stars. Such DLAs may be the most suitable environments 
for measuring the true yields from zero- or low-metallicity 
stars, free from the complications of stellar abundance mea- 
surements and the possibility that these stars are polluted 
by an unseen binary companion, or self-polluted by their 
own nucleosynthesis. The extremely metal-poor regime for 
DLAs is yet to be thoroughly explored with high-resolution 
and high signal-to-noise spectra — who knows what may lurk 
there! 
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